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Abstract 



The thermodynamic limit of certain exponential corrections to the weak cou- 
^v pling expansion of two-dimensional models is investigated. The expectation val- 

ues of operators contributing to the first order coefficient of the low-temperature 
expansion of the free energy are calculated for the order 0(e~@). They are proven 
q^ | to diverge logarithmically with the volume for non-abelian models. 



It has been recently argued that the low-temperature asymptotic expansion of two- 
Q | dimensional (2D) non-abelian models can be non- uniform in the volume Q. As an 

example, authors of ]I[ involved the superinstanton boundary conditions (BC) and have 
shown that the conventional perturbation theory (PT) produces results for invariant 
functions which differ in the thermodynamic limit (TL) from those obtained with more 
standard BC, and the difference starts at 0(1/ (3 2 ). Since the asymptotic expansion of 
invariant functions in the TL is unique, the discovery of may indicate a failure of the 
PT. However, it has been shown later that the PT with superinstanton BC diverges for 
higher orders and it has been argued that the low-temperature asymptotic expansion 
^ is non-uniform in the volume only for superinstanton BC while there are good reasons 

to believe that the PT gives true, uniform asymptotics with conventional BC @-|| 
(see, however Ref.[|]]). 

If the low-temperature expansion (which is asymptotic when the volume is fixed) is 
non-uniform in the volume, one should encounter infrared (IR) divergences (logarithmic 
in 2D) either in the coefficients of the expansion or in the remainder to the expansion. 
It has been proven that the coefficients of the continuum PT done in the dimensional 
regularization with mass regulator term is IR finite ||. In this paper we assume the 
same holds for bare lattice PT with periodic BC. This assumption however cannot 
influence our result. We are going to show that logarithmic divergences do appear in 
the exponential corrections to the finite-volume asymptotic expansion, already in the 
order 0(e~ 13 ). 

We work with £77(2) model but to make our procedure clearer we perform the 
same calculations for the XY model where it is known rigorously that the PT gives a 
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uniform asymptotic expansion |J. The arguments are very simple, both conceptually 
and technically. First, let us discuss the general strategy. In our previous papers 0, 
f| we proposed to use an invariant link formulation to perform the low-temperature 
expansion of 2D models. The starting point of calculations is the following partition 
function of 2D SU(N) model on the periodic lattice 



J\ dVi exp 



pj2 ReTrV i+ lnJ ( v ) 



where the Jacobian J(V) is a product of SU(N) delta-functions taken over all plaque- 
ttes of 2D lattice 



j(v) = n 






(2) 



The sum over r runs over all representations of SU(N), d r = Xr(I) is the dimension 
of the r-th representation. The SU(N) character Xr depends on an ordered product of 
the link matrices Vi along a plaquette 



IJ V t = V n (x)V m (x + n)V+(x + m)V+(x) . 



(3) 



The similar formulae can be written for the XY model. In either case the formal 
low-temperature expansion takes the form M 



OO -I 

fc=l ^ 



M(h,s) 



(4) 



where Aj, are known differential operators acting on the generating functional M. We 
find up to a constant |]7fl 



M{h,s) 



Y[da k (x) Y[duj k (l)exp 



"2^(0 - ivk(l)[ak{x + n) - a k (x)} 



x J^ exp 

m(x)=— oo 



2ni\J2 



mix) 



\ k 



1/2 



+ 53 u k {l)h k (l) + J2 a k {x)s k {x) 



(5) 



The group index k takes the value k = 1,2,3 for SU(2) and k = 1 for the XY model. 
The external sources h k (l) are coupled to the link fields uj k {l) while s k (x) are coupled to 
the auxiliary fields a k (x). In the abelian case s k (x) = since the operators A k depend 
only on link fields. A = y/2(3 for ST/ (2) and A = jfi for the XY model. 

The finite- volume asymptotic expansion arises from the configuration {m x } = 
for all x ("vacuum" sector). In || we have shown that the first two coefficients of the 
fixed distance correlation function calculated in this vacuum sector coincide with the 
conventional result. For general reasons (the uniqueness of the asymptotic expansion 
in the finite volume for a given BC and the IR finiteness of the coefficients of the 
expansion) one should expect that these coefficients coincide up to an arbitrary order. 



All other configurations in the sum over m x in (|5]) are exponentially suppressed with 
/3 and do not contribute to the asymptotics. We shall however prove that 



^C 1 XY = ^A 1 (d h ,d s )M(h,s) 



32 



+ 0(pe 



-TT 2 ^ 



for the XY model (actually, this is true for any A k ) while 



2I 2 C SU2^ , }[yl 



1 A 1 (d h ,d s )M(h,s) = -^ + 0({3\nLe- 2 * 2 V 



(6) 



(7) 



for 577(2) model. The constants 1/32 and 3/64 are the standard first order coefficients 
of the free energy expansion. 

Before we proceed to calculations we would like to present some simple arguments 
why Eq.(^) should be expected. In the vacuum sector the SU(N) delta-function in 
Eq.(|ID reduces to the Dirac delta-function so that the partition function (|l|) becomes 



Z(/3>1)= /n^exp 



/3^ReTrVJ 



n*c 



OJ' 



- P>k 



with the zero modes for the auxiliary fields omitted from Green functions 0. ou k (p) is 
a plaquette angle defined as 



V p = JJ v i = ex P iv k u k (p) 



(9) 



and has the following expansion 



Mp) = ^t\p) + 4^4 1} (p) + ^-A\p) + - 



On a dual lattice (see 0) the delta- function in (§) is expanded as 



(10) 



IPC** 

x.k 



da k (x) 



n ^^^- ex p 



x.k 



2tt 



^ 



X 



q=\ 1- 



x.k 
oo ,»' - X q 

x.k 



'M) 



- (2/?)«/2 



If we adopt the following definitions 



u k (l) 



d 



, oc k {x) 



d 



dh k (l) ' '"''"' ds k (x) 
the generating functional in the vacuum sector takes the form 



121 



M (h, s) = exp 



1 i 1 

-s k (x)G XjX >s k (x ) + -s k (x)Di(x)h k (l) + -h k (l)G w h k (l') 



(13) 



where Gu> and Di(x') are link Green functions (see the Appendix). These functions are 
IR finite by construction. It follows from the last equations that the IR finiteness of the 
first order coefficient requires the finiteness of the expectation value of the following 
operator 

U = lrEMx)^\x)) , (14) 

because only this operator includes IR divergent Green function G x>x >. One can show 
that the IR finiteness of this operator is equivalent to 



( E4 J W )o = 0, (15) 

where {...)o refers to the following partition function 

^o = /_"n^(Oexp[-^(/)]n5(4 0) (x)) . (16) 



x.k 



An even power of the real operator J2 x ,k ^k ( x ) calculated over positive measure can 
only be zero if the operator itself vanishes identically on all configurations contributing 
to (|i"6l). Indeed, the constraint in (]T6|) 

4 0) (x) = (17) 

can be trivially solved on the original lattice in terms of site variables 

u) k {h) = Uk{x) - ujk{x + n) , Uk(k) = uj k (x + n) - cj k (x + n + m) , 
Uk(k) = u k {x + m) - uj k (x + n + m) , u k (l 2 ) = co k (x) - uj k (x + m) . (18) 

For this solution, Y.xk^k \ x ) = 0- I 11 we have shown how Eq.(|I~5D is fulfilled in 
terms of link Green functions (see also below). 

The crucial point is that the IR finiteness of the first order coefficient requires the 
finiteness of the operator U not only in the vacuum sector but for all m x in Eq.(^|). 
For an arbitrary configuration {m x } the constraint ( |T7D becomes 

\ 1/2 

i 0) (x)) 2 j =2n\m x . (19) 

It is easy to find a non-trivial solution of (|H]) on which the operator J2 x ,k ^i ( x ) does 
not vanish (e.g., the "vacuum" solution (|18|) for the components k = 1,2 and some 
non-trivial vortex solution for ou 3 (l)). Since the integration measure is positive for any 
configuration {m x }, there is no way for Eq.(|i~5|) to be satisfied in "non- vacuum" sectors. 
This argument explains why Eq. (|7]) is not surprising but should be expected on general 
grounds. 




We want to prove now Eq.(|J). Calculating the generating functional ([5]) for the 
XY model we find (sum over all repeating indices is understood) 

M XY {h) = e L i hlG »' hl ' J2 <KE m x )e- 7r2 > 3mxG ^' m *' +7r ^ hlDl{x ' )m *' . (20) 

m, x 

Actually, Eq.(|6|) follows directly from this representation for Mxy{h): the generating 
functional depends only on functions which are IR finite. The operator A\ is given by 
j^Y^iid^/dhf). It leads to the following expansion 

A x (d h ) M X y{h) = Z vor [— - — ^ E D(x x - x 2 )(m Xl m X2 ) vor 



+ J5T2 E ^Di{xi)D l (x 2 )Di(x 3 )D l (x 4 )(m xl m X2 m xs m X4 ) vor } . (21) 

^OL/ XlX2X3X4 i 

Here, the expectation value (...) VO r refers to the vortex part of the XY model 

Zvor = E 5 (E m x )e-* 2l3m * G ^' m *' . (22) 



m x x 



The first exponential correction comes from the configuration m x = 5 xy — 5 XZ . We find 
■1 E D{ Xl - x 2 )(m xl m X2 ) vor = A £ D(y _ z)e- 2 - 2 ^y-^ [l + 0(e^)] , (23) 

L x x x 2 L y^z 

and similar expression can be written for the last term in ([H]). Since D(y — z) > 1/2 
for y ^ z, Eq.([5|) follows. 

Let us now turn to the SU(2) model. It seems to be a difficult task to calculate the 
generating functional for an arbitrary configuration {m x } because of the square root in 
(H) . We therefore restrict ourself to the first exponential correction which again arises 
from m x = 5 xy — 5 XZ , i.e. we write 

oo 

E exp[27riA E iri(x)a(x)] = 1 + 4 E cos(2tt Xa(y)) cos(27rAa(z)) + ... . 

m(x)=— oo x V z / zz 

Using the inverse Laplace transform 

1 /7/\ 1 / 2 fC+ioo A+ 

cos(2a^V /2 ) = U l ) / 4^"^ > C > ° > 

2l \1X ) JC-ioo \Jt 

the generating functional can be reduced to the form 

M su(2) (h,s)=M (h,s) (l + 4EM 1 (M+0(e- 4 ^)) , (24) 



where Mo is given in (|T3|). An integral form of Mi reads (Go = G x>x ) 

*P o -4**/3G f C+lQO dtds 



Mi(h,s) = --—e 



C-ioo (tS- l) 3 / 2 



[tGyz — Go)(sG yz — Gq) 



exp 



2ir 2 (3G yz (t + s) + 



AG yz {ts-l) 
We use the notation (hk(l) = h k (n;x)) 



(ta 2 k (z) + sa\{y) - 2a k {z)a k {y)) 



(25) 



a k{y) = E G *y [hk{x) - isk(x)} , h k (x) = J^ [h k {n; x) - h k (n; x - n)] 



n=l 



As an example, let us give an expression for Mi at vanishing sources 

Mi(0, 0) = \ e -^mv-*) (l - ^-D 2 (y -z))+0 ( e - 8 - 2 ^) . 

2 \ Gy Z J 

One sees that the convergence to the TL is very slow, like 0(1/ In L). 
The proof of Eq.(|7]) is now straightforward. The operator A\ reads |7]] 



A = ^E (e^(oV - ^E (E^(o) -«E« fc ( 

l \ k / I \ k I x k 



,( 2 ) 



where [/ is given in (14). In M we have shown that 



-L Al (d h ,d s )M (h, S ) = ^. 



(26) 



(x)^(x)-f/, (27) 



(5» 



It is seen from ( p^) that the IR divergent function G^/ can arise only from the third 
and fourth terms. It is easy to prove that the operator J2x J2k a k(x)iu k (x) is IR finite. 
Therefore, it is sufficient to study the contribution from the last term in Q2"7|). We find 



2L 2 



U {d h , d s ) M su(2 )(h, s) = U div + U fin . 



(29) 



Ufi n denotes a finite contribution which depends only on IR finite link Green functions. 
For Udiv w e find (up to the terms vanishing in the TL) 



U, 



div 



2L 2 



7 , G X:X > 



10 yz 



4 4 



where 

w x% = E E ( G w;, G ii,i;, - G u^, G i^) > 

W^Ay, z)=Y.J: (G^B v - G ia B lhV + G^B ia - G ljX B kJ >] 

i<3 v<j' 



(30) 

(31) 

(32) 
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B ltV {y,z) = (A(y) - D t (z)) (D t r(y) - D v (z)) . (33) 

The link k (E,) refers to one of four links attached to a given site x (a/). Let us 
substitute G XjX i = Go — D{x — x') in (|30D . Since the function D(x — x') is IR finite we 
conclude that the IR finiteness of Eq. (|30|) requires 



£wg, = o, £w<£(y,s) = o. (34) 

x,x' x,x' 

The first of these equations is nothing but another form of Eq. (|i~5|) . Its validity can be 
proven directly. Using (HS) we rewrite (|34|) as 



E Wi% = ^ E Mhh) [A(hb 2 ) - Aihb,)] , (35) 

x,x' b\bi 



where sums over bi run over all links of the lattice and 

A(b l b 2 ) = J2T,G bl uG b2h . 



The equality A(b\b 2 ) — A{b 2 bi) = follows from ( pip and (^5|). This equality ensures 
the IR finiteness of the first order coefficient in the vacuum sector. 

However, the second equation in (31) cannot be satisfied, in accordance with our 



general arguments. Indeed, using the same tricks as above we get 

E W®(y, z) = ~ E M - F b (z)f , (36) 

where 

F b (y) = EE( G *A,(y) " ^A,(y)) • (37) 

a; i<i 

It follows from the last two equations that in order to ensure the IR finiteness, the 
function F b (y) must be y- independent. In fact, F b (y) does depend on y as it follows from 
the definition of functions G b i and Di(y). To see this dependence directly, one can use 
Eqs.(|43|), ( f44"|) (see the Appendix) for the link Green functions and perform calculations 
in momentum space. From (|36| ) it also follows that Y, x ,x' W x x i{y, z) depends only on a 
difference (y — z). Since WQ j,(y, y) = we end up with (omitting finite terms) 

2^ E Gx,x> E e" 4 "^-*) Wj J(y, z) = Goe- 2 - 2 / 3 E <i(0, 1) + O(e^), (38) 

where c > 2. Since G ps - InL, this completes the proof of Eq. ([!]). 

We have also computed the second term in ( [30|) numerically. For /3 = 2 we find the 
values 3.036 • KT 15 for L = 200, 3.267- 10" 15 for L = 300 and 3.558 • 10" 15 for L = 500. 
These numbers are so small compared to the first order coefficient 3/128 (at (3 = 2) 
that even if there are deviations from the conventional PT in the TL, it seems there is 
no chance to see them in the Monte-Carlo simulations. 



Let us add some comments. 

I. 2D SU(N) principal chiral models are essentially compact models whose Gibbs 
measure obeys a certain periodicity required by the invariant group measure. As a 
rule, the weak coupling expansion needs a construction of the Gaussian ensemble which 
destroys the compactness and periodicity of the original measure. When the volume is 
fixed the corrections coming from the extension of the integration region to infinity and 
from the breaking of periodicity are exponentially small and thus do not contribute to 
the asymptotic expansion. The subtle question is what happens with these corrections 
in the large volume limit, in particular whether they are infrared finite or not. In 
our opinion, the link formulation of these models presents a good tool for resolving 
this problem. Summation over m x in the expression for the generating functional ([5]) 
appears due to using group invariant delta-function. In the abelian case it corresponds 
to the summation over vortex configurations. Probably, some similar interpretation 
can be given in the SU(2) case (it is, of course, also tempting to speculate that these 
configurations are responsible for a mass gap generation similar to vortices in the XY 
model). In this paper we have investigated certain exponential corrections to the first 
order coefficient in the expansion of the free energy which appear from the mentioned 
summation. Our main result, Eq.(^), shows that the IR divergences do not cancel for 
non-abelian models. 

However since such divergences appear only in exponentially suppressed terms one 
could ask how significant they are, in particular for the construction of the continuum 
limit. One considers asymptotics obtained by letting L go to infinity as a certain 
function of (3, e.g. for the XY model one takes L = 0(expv73). The proof of the 
asymptoticity then follows from certain inequalities and the fact that the PT coefficients 
with Dirichlet and free BC differ only by exponentially small corrections for sufficiently 
large f3 (see section 4 of 0). A strategy of a similar proof for non-abelian models was 
outlined in ||. Since one can take a power dependence on (3 for L, the second term on 
the right-hand side of Eq.(|7|) remains properly bounded. If however one is interested 
in the limit L — > oo where L is (3- independent, we conclude that the coefficients of the 
low-temperature expansion are not bounded uniformly in L. 

II. The second type of exponential corrections comes from the restriction of the 
integration region over link variables. For example, the one-link integral in the XY 
model is given by 



pxyP 
h = / ^du(l)exp 
J-kJb 



— uj 2 (1) — ioo(l)[a(x + n) — a(x)] 
= (2 7r ) 1 / 2 e -5(«(-+")-«(-)) 2 Lrf L^/2 + -^={a{x + n) - a{x))\ + c.c. J , (39) 

where eri(z) is the error function. From here one derives uniform bounds both for 
Mxy(h = 0), for the coefficients A k and for the correlation function |§. Our result 
actually shows that such uniform bounds in L do not exist for the coefficients of the 
asymptotic expansion of the non-abelian model, if L is /5-independent. 



Appendix 

The functions Gw and Di(x') are defined as 

Gw = 2Si t ii — G X)X i — G x+ntX i +n > + G X)X i +n i + G x+ n tX > , (40) 

Di(x') = G XjX > - G x+n . x > = D(x - x' + n) - D(x - x') , I = (x, n) , (41) 

where 

I L ~ l 1 _ e ^k n (x„-y n ) 

D(x-y) = G -G x , y = -Y: 9 v2 —WZ , ^n^O. (42) 

They satisfy the following equations 

G/^/ + Gi 2 f — Gi 3 i> — Gi 4 i' = , (43) 

D h (a/) + A 2 (a; 7 ) - D h {x') - D h {x') = 28 x>x , , (44) 

/ , (C/ 1 6 1 Gz 2 b 2 + Gi^G^b? — Gi 1 b 2 Gi 2 b 1 — Gi 3 b 2 Gi 4 b!) = , (45) 

X 

where four links U are attached to a given site x. One proves the following "orthogo- 
nality" relations for the link functions 

X; G lh G hV = 2G„, , X D b {x)G bV = , (46) 

b b 

J2Db(x)D b (x>)=2G x , x , , (47) 

6 

where the sum over b runs over all links of 2D lattice. 
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